function approx = Test2()

dt = 0.1;

j = 1;

while (dt >= 0.0001)
	i = 1;
	dx = 0.1;
	while (dx >= 0.0001)
		[p, vt, Phi, Delta, S] = PDEFD(0.2, 0.1, 50, 50, dx * 5, dt * 5, 1, 'CNS', 'put', 400, 40);
		approx(1, j, i) = dt * 5;
		approx(2, j, i) = dx * 5;
		approx(3, j, i) = p;
		[p, vt, Phi, Delta, S] = PDEFD(0.2, 0.1, 50, 50, dx, dt * 5, 1, 'CNS', 'put', 400, 40);
		approx(1, j, i + 1) = dt;
		approx(2, j, i + 1) = dx * 5;
		approx(3, j, i + 1) = p;
		i = i + 2;
		dx = dx / 10;
	end
	i = 1;
	dx = 0.1;
	while (dx >= 0.0001)
		[p, vt, Phi, Delta, S] = PDEFD(0.2, 0.1, 50, 50, dx * 5, dt, 1, 'CNS', 'put', 400, 40);
		approx(1, j + 1, i) = dt;
		approx(2, j + 1, i) = dx * 5;
		approx(3, j + 1, i) = p;
		[p, vt, Phi, Delta, S] = PDEFD(0.2, 0.1, 50, 50, dx, dt, 1, 'CNS', 'put', 400, 40);
		approx(1, j + 1, i + 1) = dt;
		approx(2, j + 1, i + 1) = dx;
		approx(3, j + 1, i + 1) = p;
		i = i + 2;
		dx = dx / 10;
	end
	dt = dt / 10;
	j = j + 2;
end

